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We revisit the renormalization of <l>-derivable approximations from a slightly different point of 
view than the one which is usually followed in previous works. We pay particular attention to the 
question of the existence of a solution to the self-consistent equation that defines the two-point 
function in the Cornwall-Jackiw-Tomboulis formalism and to the fact that some of the ultraviolet 
divergences which appear if one formally expands the solution in powers of the bare coupling do 
not always appear as divergences at the level of the solution itself. We discuss these issues using a 
particular truncation of the $ functional, namely the simplest truncation which brings non-trivial 

■ momentum and field dependence to the two-point function. 
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\Q • I. INTRODUCTION 

(-H ' Developed in the context of non-relativistic many-bodv theory [1- 3], the ^-derivable approximation scheme was 

' later generalized to relativistic quantum field theories I4| where it is considered a promising tool to address some 

O ^ of the currently open fundamental questions. This functional method relies on the two-particle irreducible (2PI) 

[ effective action, which upon a diagrammatic expansion in 2PI skeleton graphs leads to a systematically improvable 

I approximation scheme of a given quantum field theory. Remarkably, the global symmetries of the classical action are 
preserved at any order of the truncation and the solution of the self-consistent propagator equation obtained through 

^-H i: a variational principle satisfies thermodynamical consistency [2| . This feature makes the method particularly suitable 

^ ■ for the calculation of thermodynamical quantities, such as the entropy Q or the pressure [sj, which show a rather good 

, convergence behavior as the order of the approximation is increased. In an out-of-equilibrium setting, the success 

• of the 2PI evolution equations in describing the late-time dynamics of quantum fields from far-from-equilibrium 

, initial conditions {7^ has opened the possibility to study a number of problems, such as preheating ^ i9|], transport 

. . coefficients or topological defect formation [ll|, [12] . 

^ : 

. Together with the increasing number of applications of the 2PI effective action, new insight has been gained on 

■ technical aspects regarding its renormalization [T3l - [T5| . Part of the difficulty in this context comes from the fact 
. . , that the two-point function is given in terms of a self-consistent and thus non-perturbative equation, also known as 

the gap equation, which makes it difficult to identify and eliminate the ultraviolet (UV) divergences. Interestingly 
however, part of the intuition that we have gained regarding the renormalization of ^-derivable approximations 
comes from perturbation theory. In fact, one can always consider a formal (partial [l6| or complete |14 |) perturbative 
' expansion of the solution of the gap equation which allows to unveil the structure of the corresponding UV 
divergences. This structure can then be expressed in terms of certain non-perturbative objects fulfilling their own 
self-consistent equations, in particular a four-point function obeying a Bethe-Salpeter-type equation [l^ . It 
follows that the renormalization procedure for ^-derivable approximations can be formulated solely in terms of 
these non-perturbative objects, without any reference to the formal expansion used to identify the divergences, 
and it is then readily applicable in practice. However, since the previous construction is based on a perturbative 
expansion, the only rigorous statement that can be made a priori concerning the renormalization procedure is 
that, if one would formally expand the solution of the gap equation in powers of the renormalized coupling, the 
coefficients of such an expansion would converge as the ultraviolet cut-off A is taken to oo. But what about the 
solution of the gap equation itself, before it is formally expanded? Is the renormalization procedure we referred to 
above sufficient to make the solution insensitive to the cut-off when the later becomes large? Numerical studies 
seem to indicate that this is indeed the case [1, [2l|. However, proving this fact rigorously is a difficult matter 
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because the UV structure of the solution of the gap equation usually reveals itself after solving the equation and in 
general this can be achieved only numerically. Although there exist certain analytic arguments which corroborate 
these numerical observations [l3L [l^. they are all based on two important assumptions. First, that the solution of 
the gap equation exists, if not for arbitrarily large, at least for sufficiently large A. Second, that the asymptotic 
behavior of the solution is only mildly modified with respect to that in perturbation theory, that is up to some 
powers of logarithms of the momentum. These two assumptions are again difficult to prove analytically and usually 
one needs to rely on numerical evidence. A related question is how the divergences of the formal perturbative 
expansion of the solution of the gap equation appear at the level of the solution itself. Since such divergences are 
to be absorbed in a redefinition of the bare parameters, this question can also be stated as follows: what is the 
actual role played by the bare parameters in the existence and the large- A behavior of the solution of the gap equation? 

In this paper we put forward a slightly different point of view towards understanding the renormalization of $- 
derivable approximations, which might shed some light on the issues mentioned above. Generalizing to a momentum 
dependent self-energy our approach initiated in [131, discuss the behavior of the solution of the gap equation 
as the cut-off increases using a combination of numerical and semi-analytical methods, without ever relying on any 
perturbative expansion. This approach also addresses the question of the existence of a solution of the gap equation 
for arbitrary large values of A. Part of the originality of this work is that, in a sense, it revisits the question of 
renormalization of ^-derivable approximations from scratch: we shall recover known results, but from a different 
perspective. In particular, we point out a difference between the role of the bare mass and that of the bare coupling: 
if the bare mass is needed to absorb divergences of the solution of the gap equation, the role of the bare coupling is 
somewhat different, at least in the approximation that we consider here. We concentrate on the simplest ^-derivable 
approximation of the one-component ip'^ scalar field theory which introduces a non-trivial momentum dependence for 
the two-point function. The corresponding self-energy M'^{K) obeys the following gap equation: 

M^K) = ml + ^<j>' + ^ f G{Q) -^<t>' I G{Q)G{Q - K) , (1) 

|Q--fs:|<A 

where and Aq denote the bare parameters and G{Q) = 1/(Q^ -I- M'^{Q)). In what follows, the two integrals ap- 
pearing in Eq. ^ are named respectively the tadpole integral and the bubble integral. The function M'^(K) is referred 
to as the self-energy, although it includes the tree-level bare mass. It depends on the cut-off A and on the field </>, but 
we shall leave these dependencies implicit. The same remark applies to the propagator G{Q). Finally, the regulariza- 
tion that we have chosen^ is such that the norm of any momentum, including the external momentum is less than A. 

It is usually said that the self-energy M'^{K), solution of Eq. ([Ij, diverges. What is meant by this is that if one would 
let the cut-off A grow indefinitely, keeping the bare parameters fixed, the solution ]VP{K) of this bare gap equation 
would diverge. Is this really so? Note that the reason why we believe that there should be divergences comes from our 
experience with perturbation theory. More precisely, if we would formally expand the solution of Eq. ^ in powers of 
Ao, the coefficients of this expansion would diverge as A — > oo, with divergences coming both from the tadpole and the 
bubble integrals. But what about the solution M'^{K), before it is expanded? Does it diverge? And do the possible 
divergences originate from both integrals? In fact, it is not even clear that the solution exists for arbitrary large values 
of A. So, how can we even start discussing divergences? We first investigate these questions in Sec.|lTl where we argue 
that the solution of the bare gap equation exists indeed for arbitrarily large values of A and diverges. Surprisingly 
however, it exhibits a purely quadratic divergence, that is a divergence void of multiplicative or additive logarithms, 
which moreover originates exclusively from the tadpole integral. Because this quadratic divergence is independent of 
the field, we are naturally led in Sec. lIIII to absorb it in a redefinition of the bare mass TOq, without redefining the bare 
coupling Aq. We observe then that the solution of the corresponding mass-renormalized gap equation shows no more 
divergences, in contrast to what one would naively expect from the fact that the bare coupling was maintained fixed. 
Instead, for some values of the parameters, there exists a "critical" value of the cut-off above which the equation 
has no solution, and the very question of divergences is not well posed. We explain the origin of this critical cut-off 
using a "mean-field approximation" of the gap equation. This approximation is also used to illustrate another type of 
situation which could be encountered within certain ^-derivable approximations, namely the existence of a continuum 
limit which is however only reached for extremely large and practically inaccessible values of the cut-off. Both the 



^ We choose here a regularization such that the norm of any momentum is Hmited by A. This regularization sUghtly differs from the 
regularization which amounts to hmiting the norm of loop momenta only. In fact, only the former seems to be defined consistently as a 
regularization at the level of the path integral. It is also easily generalizable to smooth regulating functions, see App. 
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existence of a critical cut-off and the very slow convergence towards a possible continuum limit, although they cannot 
really be referred to as divergences of the solution of the gap equation, present the same limitation than a divergent 
solution: they prevent the definition of a cut-off insensitive self-energy. In fact, these inconvenient features can be 
traced back to the presence of logarithmically divergent contributions in the mass-renormalized gap equation. These 
divergences can be absorbed in a redefinition of the bare coupling, but we stress the fact that these are divergences 
of the equation itself not of its solution and as such, the role of coupling rcnormalization is not to absorb divergences 
of AI^{K), at least in the approximation that we consider in this work. Coupling rcnormalization is implemented 
in Sec. IIVI and the corresponding completely renormalized gap equation exhibits a solution which seems to exist for 
arbitrarily large values of the cut-off A, at least in some relevant range of parameters, and converges towards a certain 
limit as A ^ oo. This limit is approached to a very good accuracy (oc 1/A when using a sharp cut-off) already for 
"accessible" values of the cut-off. 



II. BARE GAP EQUATION 



Let us first assume that the bare parameters uiq and Aq are kept independent of A. We refer to the corresponding 
gap equation as the bare gap equation. Our goal is to discuss the behavior of its solution M'^{K) as A is increased 
for fixed K. To this purpose, it is convenient to consider a more general situation where some of the components of 
the 4-vector K grow with A as well, that is K = AK + L, with \AK + L\ < A. Performing the change of variables 
Q = AQ and introducing the notations M^{Q) = AP{Q)/A^, mg = m§/A^ and (jP' = (j)^ / A^ , the bare gap equation 
^ becomes 



Ao 
2 



IQKi 



G{Q)G{Q -k- L/A) 



(2) 



\Q- 



<i 

K-L/A\<1 



with G(Q) = ll{Q^- + ]VP{Q)). 
asymptotic behavior 



Because TOq 



and 



as A — )• oo, this equation is "compatible" with the 



M'^(k + L/A) ^ Af^ > 0, 



as A — )■ oo, for fixed K and L, where fulfills the self-consistent equation^ 



= — 



\Q\<i Q"^ 



M2 

-'"oo 
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(4) 



By "compatible" we mean that if one assumes that the asymptotic behavior ([3]) is obeyed by the self-energies in 
the right- hand-side of Eq. ([2| , the later produces exactly the same asymptotic behavior for the left-hand-side of the 
equation. In particular, we obtain M'^{K) = M'^{K / A)A^ ^ A^ that is, the solution of the bare gap equation 
exhibits a pure quadratic divergence at leading order in the asymptotic expansion as A oo for fixed K. Note that 
we have not proven that Eq. ([2|) admits a solution for sufficiently large A. However, the fact that the contribution 
of the tadpole integral dominates over the contribution of the bubble integral makes this plausible. This fact is 
confirmed numerically for the smallest momentum stored on the grid \K\ — km — 5.10^"* (we also take mg = 0.01 and 
(f)^ = 0.1) together with the announced asymptotic behavior, see the left panel of Fig. [TJ In order to save computer 
time, the parameters were chosen such that the asymptotic behavior ^ is observed for reasonable, that is to say, 
not very large values of the cut-off. 



It is interesting to compute the first correction to the asymptotic behavior M^{K) ~ Af^A^. This is done in 
App. m We obtain: 



Ao 



X2 

,2 



GliQ) 



X X2 



QI<1 



GliP) / Goo{Q)Goo{Q - P) . 

P\<l J\Q\<1 
\Q-P\<i 



^ If we consider each side of this equation as a function of M^, the left-hand-side increases hnearly from to oo, whereas the right- 
hand-side decreases from a strictly positive value to 0, if Aq > (which we assume to hold throughout this section and the next one). 
Then, because both sides of the equation depend continuously on , there is always a unique and strictly positive solution , when 
Ao > 0. 
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1: The left plot represents the ratio M {K)/A as a function of logjg A (and as a function of A in the inset) for 
= km = 5.10"* (the smallest momentum stored on the grid) and the comparison to the asymptotic estimate M^, determined 



from Eq. Q. The right plot represents the difference M^(fem) — A^ as a function of logjg A (and as a function of A in 
the inset) and the comparison to the asymptotic estimate, as given by the right-hand-side of Eq. ([S]). The parameters are 
mo = 0.01, Ao — 1, and (f)^ = 0.1. For a given A the number of discretization points is [lOOA]. 



(5) 

with Goo(Q) EE 1/{Q^ + M^) and 1/Aoo = 1/Ao + (1/2) jj^i^.^ G^((5). Equation ^ shows that there are no 

logarithmic divergences in M'^{K) as A — > oo for fixed K. This result is confirmed numerically, see the right panel 
of Fig. [TJ Note also that the divergence of M'^{K) comes from the tadpole integral only. In contrast, if we were 
to consider the perturbative expansion of M^{K) in powers of Aq, the corresponding coefficients would diverge with 
divergences originating both from the tadpole integral and from the bubble integrals. This result illustrates that, 
upon resummation, some of the divergences which appear in the formal perturbative expansion of the solution of the 
gap equation, do not appear in the solution itself. Does this mean that there is no need to absorb these perturbative 
divergences? This is one of the questions we shall investigate in the next sections. 



III. MASS-RENORMALIZED GAP EQUATION 

The quadratic divergence M^A^ discussed in the previous section does not depend on the field cf) and we can then try 
to "absorb" it by adjusting the bare mass TOq, which then becomes a function of A and keeping the coupling Aq fixed. 
As already mentioned, the divergence originates from the tadpole integral. However, because the latter depends on 
the field 4> through the solution M^(Q), it cannot be completely absorbed in a redefinition of TOq. Still, we shall see 
below that, in order to remove the divergence, it is enough to choose ttiq = ^(Ao/2)[/|q|^^ Go((5) + finite part] with 

Go{Q) = + 171^) and some fixed renormalized mass. The bare mass is defined up to an additive "finite part" 

which we can choose such that the renormalization condition = at (jP = is fulfilled."* This amounts here to 
the particular choice: 

ml = m'-^ f GoiQ) . (6) 

^ "'|Q|<A 



^ Due to the presence of the term Ao0^/2 in Eq. ((TJ, renormalizing the coupling Ao would mean that the divergence of M'^{K) contains 

a terms proportionnal to ijp , which it has not, see Eq. (|3]l 
4 For 4? = 0, the self- energy becomes momentum independent and we do not need to specify a particular momentum in the renormalization 

condition. Other choices of renormalization points are possible. In particular, if one is interested in discussing the broken phase, 

one should rather impose a renormalization condition at a given non-vanishing value of the field, or at vanishing field but non-zero 

temperature. 
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Plugging this expression into the bare gap equation, we end up with 



2 I •^0,2 



Ao 
2 



IQKA 



G(Q)-Go(Q) 



'^0 i2 



G{Q)G{Q~K). 



(7) 



'IQKA 
|q--r:|<a 



which we refer to as the mass-renormalized gap equation. Note that, because we have modified the dependence of 
the equation with respect to A, we need to analyze once again the behavior of the solution of Eq. ([7]) as A increases. 
Using the same rescaling and notations as in the previous section, we obtain 



M\K + L/K) ^ + 



+ 



An 



|Q|<1 



G{Q)-Go{Q) 



'^O 72 



G{Q)G{Q -K- L/A) , (8) 

|Q|<1 
K~L/A\<1 



where TO^ = ^s/A^ andGo(Q) = 1 / {Q^ + {Q j) . This equation is compatible with Af2(X+ L/A) = M^{AK+L)/A^ 
becoming smaller and smaller as A increases for fixed K and L (this will be confirmed numerically in what follows) and 
supports then the fact that the choice © for the bare mass eliminates the quadratic divergence. Notice however that 
the question of the existence of a solution for any value of A is now more subtle than in the case of the bare gap equation 
because the contribution of the subtracted tadpole integral does not necessarily dominate over the contribution of 
the bubble integral. We shall now study more in detail what happens to the solution of the mass-renormalized gap 
equation as A is increased. 



A. Effect of the bubble integral 

Before discussing the mass-renormalized equation ([T]), let us consider a simpler equation where we focus on the 
effect of the bubble integral only, that is: 

M\K) = ml-^f G{Q)G{Q~K), (9) 

^ J\Q\<A 
\Q~K\<A 

where = + (Ao/2) 0^ and Aq(/)^ are considered as independent parameters. Such a simplification makes possibly 
sense when the field (f> dominates over the mass m. Its main purpose here is to introduce, in a simpler context, the 
arguments and tools that will be used in order to deal with Eq. ([7]). Moreover, some of the properties of the solution 
of Eq. ([7]) will be visible here already. 



1. Numerical solution 



For reasons that will become clear below, we solve Eq. © using a differential method, that is by increasing the 
cut-off A continuously from the initial value A = at which the solution of Eq. ([9]) is known and given by rrv^. An 
"infinitesimal" change from A to A -I- SA is encoded in an evolution or flow equation which is more conveniently derived 
for M^{AK) at K fixed. ^ In fact, after a simple change of variables, we can write 



M^{AK) = mi 



\Q\<1 

)-k\<i 



G{Q)G{Q-K), 



(10) 



where we note that the integration domain does not involve the cut-off A anymore. Then, taking a derivative with 
respect to A at fixed K , we arrive at the equation 



5aM"(AX) = \l<t>'^ 



G\Q)G{Q~K) 



Q\<1 

k\<i 



±dAM\AQ) - ^M'(AQ) 



(11) 



^ If A 7^ 0, it is completely equivalent to consider R'P{K) for \K\ < A or RP{AK) for \K\ < 1. 
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which can be solved from the initial condition^ liniA^o M'^{KK) 



using a Runge-Kutta algorithm combined 



with a linear solver based on LU decomposition (see App. |E] for details). Note that Eq. pT|) is compatible 
with limA->o d!yM'^{KK) = 0. The flow is then flat at initialization. For practical purposes it is then more 
convenient to choose a non-zero, but small value for the initial A. We have chosen Ajnit — 2krm where 
is the value of the smallest momentum stored on the grid, which is kept fixed in a given run. The result of 
integrating the flow equation is shown in Fig. [2] for a particular choice of the parameters and Agc/)^. Each curve 

corresponds to M'^{K) for a given value of the cut-off A, plotted as a function of the rescaled momentum |A'| = |A'|/A. 
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FIG. 2: Evolution with A of the solution M'^{K) of Eq. (fTTI) plotted as a function of the rescaled momentum \K\ = \K\/A. 
Lower curves correspond to increasing values of A. The flow cannot be continued above a certain critical cut-off Ac. The plain 
curve represents the solution obtained as A — )■ A~. The dashed curve represents the limiting function p2|) . which is approached 
by M'^{AK), at least for \K\ close to 1 and if Ac is large enough, see the text for details. The parameters are — 1 and 
\o4>^ — 50, the smallest momentum stored on the grid is km ~ 7.5 10~* and the number of discretization points is 175. 



We first observe that as A increases the self-energy (A AT) becomes insensitive to the cut-off in an increasing range 
of I AT I, starting from the highest value \K\ = 1. This can be understood easily from Eq. ([TU|) . which is compatible 
with M^{Q) becoming smaller and smaller as A increases and thus with A/^(AAr) approaching 



1 



1 



QKi QHQ 

'K\<1 



Kf 



(12) 



This is confirmed in Fig. [51 Of course, since the integral in diverges logarithmically as AT — >■ 0, we expect the 
agreement between M^(AAr) and (|12l) to be the best for \K\ close to 1 and to extend over a wider range of |Ar| as A 
is increased, which is also visible in Fig. [2] As we shall see in a moment, there exists an upper value of the cut-off 
which cannot be overpassed. This explains why in Fig. [21 where this upper value is only slightly larger than 12.5 
(m^ — 1 and Agc/)^ = 50), the agreement between M^(AAr) and (|12p is only effective in the region where \K\ is really 
close to 1. 



The previous remarks do not apply to the neighborhood of AT = where Af^(AAr) remains sensitive to the cut-off 
A. In order to study more precisely what happens in this small momentum region, in the left plot of Fig. [3| we consider 
the evolution with A of the self-energy evaluated for \K\ = 0,*" which we denote by = M'^(O) from now on. We 



® This limit is obtained from Eq. I|10|l for fixed K using a consistency argument. 
For numerical convenience we rather evaluate M^(K) at \K\ = fcm, which should be regarded as our numerical approximation of J\/^(0). 
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observe that, after a relatively smooth evolution, the derivative of with respect to A becomes infinite (numerically 
at least) for some "critical" value Ac and the mass approaches a non zero value ^0. A fit in the vicinity of 
Ac reveals that the singularity is well approximated by 

M^-M^oc (Ac-A)'^'. (13) 

The flow equation suggests that this singular behavior "propagates" to non-zero values of K, as we confirm in 
the right plot of Fig. [3] for the highest rescaled momentum available on our grid, that is \K\ — 1. Note also that, 
despite the fact that the smallest value of the self-energy, obtained at zero momentum, decreases with increasing A 
(due to the bubble integral contribution), the flow stops before the propagator develops a pole at zero- momentum. 




I ^ , , , ^ ^ 1 0.91 I ' ' ' ' ' ^ 1 

2 4 6 8 10 12 14 2 4 6 8 10 12 14 

A A 
FIG. 3: Left plot: evolution with A of the self-energy M'^{AK) evaluated for \K\ = (numerically \K\ = fcm/A). Right plot: 
evolution with A of the self-energy M^{K) evaluated for the highest \K\ = 1. The parameters are — 1 and Ag<;/>2 = 50, 
while fcm — 7.5 10"'*, and the number of discretization points is 175 for the solution of the flow and [400A/3] for the iterative 
solution. 



The presence of a singularity in the flow does not necessarily mean that there is no solution to Eq. ^ above Ac . It 
only means that if a solution exists for such values of A, it cannot be accessed by integrating the flow equation from 
A = 0, without any additional information. For a given A, we can also solve Eq. Q using iterations from an initial 
(constant or perturbative) ansatz for the self-energy or from a solution to the gap equation obtained at a different 
value of the cut-off. The results of such an iterative procedure are represented in the plots of Fig. [3] for some values of 
A and they compare pretty well with those obtained from the flow equation. We observe that the iterative procedure 
fails around Ac as well. This again does not necessarily mean that there is no solution to Eq. ^ above Ac but only that 
if a solution exists, it cannot be reached within the iterative approach from the initial ansatz that we have considered. 
In order to better grasp the origin of the singularity, we shall now consider a kind of "mean-field approximation" 
for Eq. ^ which, because it will decouple the different momenta, will turn the original integral equation ^ into 
a certain number of simple numeric equations. The mean-field approximation will not only allow us to understand 
analytically the origin of the singularity, it will also allow us to argue that there is indeed no solution to Eq. © above 
Ac. A similar approximation will be used later to discuss the properties of the original mass-renormalized equation ([?])■ 

Before proceeding further with the construction of a mean-field approximation, we add a few remarks concerning 
the numerical implementation. On the one hand, since the integration of the flow equation ([TT|) involves the (time 
consuming) numerical resolution of a linear system, we are much more limited concerning the number of discretization 
points in the case of the flow approach, than in the iterative one. Consequently, the latter provides a more accurate 
solution, even though discrepancies with the flow approach only start to appear very close to Ac. On the other hand, 
the iterative method is not the best choice for approaching Ac, since the slightest overpassing of Ac will only be visible 
after an important number of iterations have been considered and leads thus to an important slowing down in the 
process of finding Ac . What is actually observed, by monitoring the change of (0) between two consecutive iterative 
steps, is that the iterative method seemingly converges in its first stage, but after a large number of iterations the 
procedure clearly diverges. By increasing the precision of the numerical integration and the number of discretization 
points we could convince ourselves that it is unlikely that this feature is the consequence of error accumulation. In 
contrast, the integration of the flow equation with adaptive step-size is such that the flow remains always below Ac. 
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There is also a numerical slowing down as we approach the singularity, but it is globally less important than with the 
iterative approach, because none of the generated data is wasted. 



2. Mean-field approximation 



Let us first consider the small \K\ region. As we already argued and as we checked numerically, M'^{Q) becomes 
smaller and smaller as the cut-off increases. Then, the integral in Eq. pU)) in the region of small \K\ is dominated by 
small values of |Q| and it makes sense to consider the following approximation for the zero momentum self-energy 

^ M\0) ^ml-^l , (14) 



IQKA 



(Q2 + Af 2)2 



obtained from Eq. ^ by setting K = Q and replacing the self-energy appearing in the bubble integral by its value at 
Q = 0. For non-zero values of \K\, we consider the following approximation 

Q-A'|<A 

where we have replaced the integral by the contributions obtained when one of the propagators reaches the smallest 
mass, that is NP = M'^{Q).^ There are two such contributions, one corresponding to Q = the other to Q + K = Q 
which should be added if we were really restricting the integration domain. But because we continue integrating over 
the whole domain, the two contributions need to be averaged. However, due to the particular UV regularization that 
we have considered, the two contributions over which we average are identical and we finally end up with Eq. (jlSp . 
As announced, we have now a certain number of decoupled equations parametrized by K. All equations depend on 
M^, which needs to be determined first from Eq. (|14l) . This approximation can certainly not allow us to access the 
actual solution of Eq. ^ and the actual value of Ac. We are however only interested in the properties of M^{K) 
as A increases and those seem to be pretty well captured by the approximation, as we can see in the comparison 
depicted in Fig. |31 Although this is not directly visible on all curves, they all present a singularity of the type 
for some critical cut-off Ac. For identical values of the parameters, the value of the critical cut-off is different between 
the "exact" equation and its mean field approximation. However, note that the corresponding evolutions with A are 
rather close to each other when plotted against the reduced cut-off A/ Ac. 




FIG. 4: A comparison between the mean-field (lines) and exact (points) solution of the respective equations (|9]) and (|14|) - (|15|) 
for M^{KK) as a function of the reduced cut-off A/Ac, at \K\ = A and \K\ — (inset). For parameters m| = 1 and Ao<^^ = 50 
the critical cut-off is A^ = 12.54 in the exact case and Ac = 9.35 in the mean-field case. 



The self-energy appears to be a monotonously increasing function of \K\, see also the discussion in App. iDl 
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Let us now use the mean- field approximation to understand analytically the behavior of M'^{KK) as A increases. 
We shall focus on Eq. ([H)) , since the appearance of a singularity in the flow of the zero-momentum self-energy triggers 
the appearance of a singularity in the flow of M'^{hK) for any other K. Our goal is to discuss how the solutions of 
Eq. ([T3]) evolve with A. To this purpose, we rewrite this equation as = fiy{M'^) with 



/a(M') 



(16) 



and study the shape and the zeros of fA{M'^) as A increases. For A = 0, fo{A'P) — ~KP -I- and Eq. (fT4|) has 
one solution only: M"^ ~ m|. To treat the case A > 0, we consider the first and second derivatives of fA{M'^) with 
respect to JVP: 



f;{M') = -3a; 



2 j,2 



(17) 
(18) 



Because fj[{M'^) < , f'i^{NP) decreases strictly from /a(0+) = oo (the integral in Eq. (fTT]) is infrared divergent) 
to /a(oo) = —1 (the same integral is suppressed for ^ A'^). If we denote by M'^{K) the only value of M"^ at 
which /^(M^) vanishes, we conclude that J^iM"^) has a maximum for = (A). More precisely, it increases 
strictly from /a(0"'') = —oo (the integral in Eq. (fT6|) is infrared divergent) to fA{M^{A)) and then decreases towards 
/a(cx)) = — oo (the same integral is suppressed for AP :s> A^). It follows that the number of solutions of the equation 
= fAiM^) depends on the sign of /A(Me (^))- To discuss this sign, note that by definition = /^(Af^ (A)). It 
follows then that 



_d_ 

dA 



M|(A) 



1 



A3 



/;'(M2(A)) dA Sn^fjliApiA)) (A^ + M2(A))3 



>0, 



and 



dA 



fA{MUA)) = 



A3 



dfA_ 

OA 167r2 (A2 + M2(A))^ 



< 0. 



(19) 



(20) 



Thus M2(A) increases strictly with A whereas /a(M2(A)) decreases strictly with A. Let us now determine the 
extremal values of ^^'^ fA{M^{A)). After performing the integral in Eq. (fT7|) . the equation defining Mg (A) 

reads 



A202 







A^ 



327r2 Af2(A)(A2 + M2(A))2 



Now, because Af2(A) increases strictly with A, it has a limit as A — > 0+, and it has a limit or goes to cx) as A 
From Eq. (PT|) . it is easily seen that the limit as A — >• 0+ is necessarily and that Af2(A) cannot diverge as A 
We finally obtain 



M^{A) ^0 as A ^ 0+ and M^{A) 

We need to be more accurate concerning the behavior of (A) as A 

Af2(A) 



K<f> 



327r2 



as A — oo . 



0+. It is easily checked that 



A2 



-J' oo as A 0^ 



Similarly, fA{M^{A)) reads explicitly 



/A(A?2(A))^-A?2(A) + m2 



3£ 

327r2 



In 



A2-hAf2(A) 



A2 



Af2(A) 



Using Eqs. (|22|) and ([23l) . we finally obtain that 

fAiM^iA))^m'' as A ^ 0+ and /A(Af2(A)) 



327r2 



A2 + Af2(A) 



In A as A -)■ oo . 



(21) 

■> oo. 
■> oo. 

(22) 
(23) 
(24) 

(25) 
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To summarize /a(M2(A)) decreases strictly from fQ+{M^{0+)) = > to f^{M^ (oo)) = -oo. Thus there 
exists a critical value Ac such that f\^{M^{Ac)) = and above which = /a(-M^) has no solution, as announced 
earlier. Notice that what drives f\{M^{A)) to negative values is the presence of a logarithmic divergence. This is 
not a logarithmic divergence of the solution (the later ceases to exist at Ac) but a divergence of the equation itself. 

The behavior around Ac can now be understood from the following argument. The critical cut-off Ac is defined as 
the particular value of A at which the maximal value /a(Mo (A)) equals zero. Denoting by the value of M^{A) 
at A = Ac,^ we have = f'/^^{M^) and = fjy^{M^). If we now expand the equation = fi^{NP) around A = Ac and 
Af2 ^ ]^,^2^ obtain 

= /a(M2) ^ /^^(Af2) ^ fl(M^)(M^ - m2) + lfl(M^){M^ - M^f + . . . 



dA 



(A - Ac) + ... , (26) 

Ac,M? 



where in the first line we have expanded up to second order because the first two terms vanish. We then obtain that, 
in the vicinity of the critical point 



I ^1 [ri{M^)dA 



(Ac - A)i/2 , (27) 



which is similar to the behavior given in Eq. (|13p . 

It is finally interesting to study how the critical cut-off depends on the parameters. It is convenient to introduce 
the rescaled quantities (j) = •^o0/(47r\/2), rh = m^/ip, Mc = Mc/(f>, and Ac = Ac/0, in terms of which the equations 
defining the critical point, that is = /j^ (^'^c) ^^'^ ^ /Ac(^'^c)j can be rewritten in the form 

= - Af2 + Afc + ln(l - Mc) , Ac = ^i^— - . (28) 

(1 - Afc) ' 

A straightforward analysis then leads to 

M^_l^_e-™' and Ac ~ e™'/^^ (29) 

as — > oo. The complete dependence of Afc and Ac with respect to rh^ as obtained from Eq. ([25)1 is plotted in 
Fig. [5] and compared to the asymptotic estimates given in Eq. ([29)) . We see that the asymptotic forms given in 
Eq. ([2^ are pretty good approximations of the exact solutions of Eq. even for small values of rh. 




FIG. 5: The dependence on the rescaled parameter rh of the rescaled critical mass Mc and the rescaled critical cut-off Ac (inset) 
as obtained from Eq. (|28p in comparison to the corresponding asymptotic estimates (|29p . 



This is also the value of the (unique) solution at A = Ac. 
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B. Combined effect of the tadpole and the bubble integral 

We now move to the original mass-renormalized gap equation given by Eq. ([7]) or Eq. ([8]), which includes both 
the effect of the tadpole and the bubble integrals. We first solve this equation for increasing values of A using a 
flow equation that we derive along the same lines as in the previous section, see the text around Eqs. pUj) and (|lll) . 
Starting from Eq. ([7]) one obtains 



Ql<i 



'IQKi 

\Q~k\<i 



&{Q)G{Q-K) 



A2 



<i 



(30) 



where GoiQ) = + and — m^/A^. This equation is solved with the method outlined in App. lEland its 

solution is presented in Fig. [51 We observe that the behavior for increasing A is similar to the one obtained previously 




FIG. 6: Evolution with A of the solution M'^{K) of Eq. H30I) plotted as a function of the rescaled momentum \K\ = |Jsr|/A. 
The parameters are = 0.49, Ao = 20, and 0^ = 4, while fcm = 3.75 10~* and the number of discretization points is 500. 



without the inclusion of the tadpole integral and that there is again a critical value of the cut-off which cannot be 
overpassed. In order to save computer time, the particular values of the parameters given in the figure caption were 
chosen such as to observe this behavior for not too large values of A. A similar behavior is observed for smaller values 
of the coupling, but then the value of Ac is larger. An analysis of the dependence of the critical cut-off with respect 
to the parameters is given below. The region around \K\ = 1 is well approximated by (jl2p even tough there is now an 
additional (momentum- independent) contribution originating from the tadpole integral. The region around \K\ — 
is still sensitive to A and, as shown in Fig. [3 close to the critical value of the cut-off Ac at which the flow becomes 
singular, the solution of the flow is approximated by Eq. (jl3p . As before, this singularity propagates to any value of 
K, and again we can try to understand the origin of the singularity using a mean-field approximation, which consists 
in an equation for the zero momentum self-energy 



Ao 
2 



IQKA 



1 



1 



Q2 + Af 2 Q2 



X2 
^0 



1 



as well as a set of decoupled equations for each value K (the equations are coupled to though): 



Q|<A 



"^0 .2 



1 



IQKA Q^ + m{Q + KY + m{K) ■ 

+K\<A 



(31) 



(32) 
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Again, we do not expect these equations to reproduce the actual values of M'^{K) and Ac. However, they seem to 
describe the large A behavior correctly, at least for small values of \K\, see Fig.|Sl Although this is not always visible, 
all evolutions plotted in Fig. [8] present a singularity for some critical cut-off. Both in the case of the exact equation 
and its mean field approximation, the singularity is approached from above for small values of |i^| and from below 
for values of |/ir| close to 1. Note also that for \K\ <^ 1, the agreement between the exact equation and its mean field 
approximation is remarkably good in the whole cut-off range, once the cut-off has been rescaled by the corresponding 
Ac. The quantitative discrepancies observed for non-small values of jii'l and A close to Ac can be traced back to the 
fact that replacing the self-energy by the zero momentum self-energy in the subtracted tadpole integral is not such a 
good approximation as in the case of the bubble integral. Below, we will discuss another limitation of the mean-field 
approximation. 




FIG. 7: Left plot: evolution with A of the self-energy M'^{AK) evaluated for \K\ = (numerically \K\ = km/A). Right plot: 
evolution with A of the self-energy M'^{K) evaluated for \K\ = 1. The parameters are = 0.49, Ao = 20, and = 4, while 
km = 7.5 10"*, and the number of discretization points is 400 for the solution of the flow and [lOOA] for the iterative solution. 




FIG. 8: A comparison between the mean-field (lines) and exact (points) solution of the respective equations for M^{AK) as 
a function of the reduced cut-off A/Ac, at various values of \K\. For parameters = 0.49, Aq = 20, and 0^ = 4, the critical 
cut-off is Ac — 89.01 in the exact and Ac — 102.73 in the mean-field case. 
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Equation (pT|) can be put in the form = g\{M'^) with 

It can be studied analytically along the same lines as the equation — /a(-M^) above. The shape of g\{M^) is 
similar to that of /a(M^): it increases from gA{0^) — —oo to (7a(M^(A)) and then decreases to —oo. The situation 
seems thus pretty similar to the one concerning the function f\(M'^): the existence of solutions to the equation 
= (7a(M^) depends on the sign of gA{M^{A)). However, the discussion of the sign of gA{M^{A)) is a little bit more 
involved than that of f\{M'^[h.)) and for this reason we relegate it to App. [C] and focus here on the results. The 
main difference with respect to the case of the function /a(M^) is that the behavior of the solution at large A 
depends now on the choice of parameters: 



In 



A^ 



r, In ■ 



M' 



K4> 

327r2 



In- 



A " + 



A^ 

A2 + A/2 



(33) 



1. If m? < Ao0^, there is a critical cut-off Ac above which the mean- field approximated equation pil) has no 
solution. This situation is similar to the one we have discussed in the previous section. The behavior in the 
vicinity of Ac is again — cx (Ac — A)^/^ which is also what is obtained by fitting the solution of the 
original equation, see Fig. [T] This comes as no surprise since the explanation for the critical behavior is exactly 
the same as in the previous section, see App. [C] The plot of Fig. |9] gives an idea of how Ac depends on the 
parameters. The critical cut-off increases rapidly as the coupling decreases, so one could argue that for moderate 
values of the coupling it plays no role. However, as shown in Fig. (TU] there is an important variation of in 
the range [0, Ac] which prevents the definition of a cut-off insensitive solution. One could try to obtain a less 
sensitive solution by restricting the interval over which A is varied. For instance, if we suppose that m — 100 MeV, 
(fy^ /m? = 0.1, and Aq = 12 then Ac/rn ^ lO**", whereas the highest relevant scale we can conceive so far is 
the Plank mass, which would correspond to A/m ^ 10^''. Even if we would restrict the cut-off to vary in this 
window, we would observe an important variation of the mass, which prevents the definition of a cut-off independent 
result. If we consider as an upper cut-off the electroweak scale, A/m ^ 10^, the variation is smaller but still important. 




5 10 15 20 25 30 
Ao 



FIG. 9: Iso-(Ac/m) curves in the (Ao, Ao0^/?n^) plane. The label on the curves indicates the value of logj^g Ac. 



2. If m? > Aq^^, Eq. (I5T]) has always two-solutions. The one which is continuously connected by the flow to the 
unique solution at A = 0, and which thus corresponds to the solution of the niass-renormalized gap equation plotted 
in Fig. |6l converges to the limit — Aq^^ as A ^ oo. As we show in App.O this limit is approached very slowly: 

InA^ \ J 
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We also show that in this case, the solution M'^{K) of Eq. (l32l) is defined for arbitrary large values of A and that the 
equation is compatible with M'^{K) having a continuum limit M^{K), solution of the following finite equation 



)l2 

-^0 ,2 



1 



Q 



1 



(Q-Kr + MliK) 



(35) 



For all practical purposes, the very slow convergence towards the continuum limit prevents the definition of cut-off 
independent result. In fact, this continuum limit cannot be considered too seriously because it is only reached for 
cut-off scales far beyond the regime of applicability of the model. As already mentioned, one could then try to define 
cut-off independent results by restricting the range of variation of the cut-off. But as it was already the case in the 
previous point, the variation of with respect to the cut-off remains important. 



2.0 n — • — • — • — • — I — • — • — • — • — I — • — • — • — • — I — • — • — • — • — r 




FIG. 10: Solution of the mean field approximated equation (|3ip for — 1, (f)^ = 0.1 and different values of the coupling 
Aq. The figure illustrates the transition between the cases 1. and 2. discussed in the text. Starting from a situation where 
< Ao'^'^, and thus such that the equation has no solution above a certain Ac, the coupling is lowered down to a situation 
where > Xo4>^, and thus such that a continuum limit exists. As we decrease Ao and approach the case = Ao<^', Ac 
increases. Right when — Xo(f>^, Ac is still finite but incredibly large. As soon as > Ao<^^, there is no critical cut-ofF 
but rather a continuum limit (thin horizontal line) which is approached very slowly, as indicated by the dashed line which 
represents the asymptotic estimate given in Eq. (|34|l . 



In contrast to case 1., it is difficult to test numerically whether case 2. is relevant for the original mass renormalized 
equation ([7]) simply because we cannot access such incredibly large values of the cut-off where a continuum limit 
could be observed. Even though case 2. could occur in certain ^-derivable approximations (see jl7| for an example 
of truncation where it occurs) , we do not believe that it occurs in the original mass-renormalized gap equation of the 
present ^-derivable approximation, as we now explain. Note first that the existence of a continuum limit in the mean 
field approximated equation can be understood as follows. Using similar manipulations as in App. |Bl Eq. pip can be 
rewritten as 



A 



"^0 ,2 



A 



IQI<A 



iM'-m'rGt,iQ)GiQ) 



AqA ,2 



IQKA 



QI<A 



GliQ) 



G\P).. 



(36) 



P|<A 



with 1/A = 1/Ao + (1/2) /|Q|<A^o(Q)- I* easy to check that this equation is compatible with the existence of a 
continuum limit for Af^. Indeed, using 



<A 



167r2 



convergent , 



(37) 
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and the fact that A ^ 327r^/ In A^, one checks that the divergences in the last two contributions of Eq. (p6)) compensate 
and that the only possible continuum limit is — Ao0^ when > XocfP, as obtained in App.O The cancellation of 
divergences is due to the fact that the last integral in Eq. (|36| is the product of two decoupled single bubble integrals. 
It follows that its leading divergence is the square of the divergence of the bubble integral: 

/ GliQ) [ G'iP) = —^(\n^)L^)+.... (38) 
J\Q\<A J\P\<A (167r2)2 V m^J \ M^J 

If we now use the same approach to discuss the original mass renormalizcd equation ([7]), we obtain 



^ -'IQKA 



- f G{Q)G{Q^K) + ^^' f GliQ) f G{P)G{P-Q), (39) 

^ J\Q\<A 4 J|Q|<A J|P|<A 

|Q--R:|<A \P-Q\<A 

where we observe in particular that the momenta in the last integral are now coupled. This coupling is enough to 
modify the leading divergence of the integral^" which is not anymore the square of the divergence of the bubble 
integral but only half of it: 

GUQ) [ G{P)G{P-Q) = l—^(\n^) + convergent, (40) 
IQKA ^|P|<A 2 [IGn^y V / 

\P~Q\<A 



as we show in App. |Dj The divergences of the r.h.s. of Eq. (p9l) do not cancel anymore and the right-hand-side 
behaves as — (Ao0^/(327r2)2)(ln A^/m^)^. This is incompatible with the existence of a continuum limit since 
the solution of the gap equation cannot become negative. It is then most probable that the solution of the 
original mass-renormalized gap equation ceases to exist beyond some value of the cut-off, just as described in 
case 1. Anyway, irrespectively of the fact that case 2. is relevant for the original mass renormalized equation 
or not, it is clear from our analysis that renormalizing the mass without renormalizing the coupling, although it 
does not lead to a divergence of the solution of the mass equation, is not sufficient to ensure the insensitivity of 
the solution with respect to the cut-off. In the next section, we cure this problem by means of coupling renormalization. 



IV. COMPLETELY RENORMALIZED EQUATION 



The analysis of the previous section shows that the presence of logarithmic divergences in the mass-renormalized 
gap equation does not always translate into logarithmic divergences of its solution. Instead, these divergences are 
responsible for the appearance of a critical cut-off A^. above which the gap equation has no more solutions or could 
also lead to the existence of a continuum limit which is however only reached for incredibly large, and practically 
inaccessible, values of the cut-off. These features prevent the existence of solutions of the gap equation which 
are insensitive to the cut-off. In this section, we explain how to get rid of these effects and define a completely 
renormalized gap equation whose results are insensitive to the cut-off already for moderate values of the cut-off (up 
to terms of order 1/A when one uses a sharp cut-off). The idea is to use coupling renormalization, not to absorb 
subdivergences of the self-energy M'^{K) (there are no such divergences as we have seen), but to eliminate the 
remaining divergences of the gap equation. In this way, one can hope to get rid of the undesirable effects mentioned 
above. 



We assume that the self-energy does not play a role in the divergences of the bubble and the double integrals. This comes from the fact 
that, in the present truncation, the self-energy does not modify the leading UV contribution to the inverse propagator and also from 
the fact that the bubble integral has a vanishing superficial degree of divergence and that the only divergences in the double integral 
come from \P\ — > oo or \P\, \Q\ oo, but not from |Q| oo. 
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A. Renormalization method 



The main difficulty is tliat one sliould be able to perform the renormalization program outlined above using a field- 
independent bare coupling. It is not necessarily possible to do so if one keeps the original form of the gap equation. 
We know how to proceed if we slightly modify the original equation into 



Ao 
2 



IQKA 



GiQ)-Go{Q) 



A2 
2 



GiQ)G{Q~K)., 



(41) 



'IQKA 

\Q~K\<A 



by allowing the bare couplings Aq and A2 to be different. Note also that the bubble integral is multiplied by 
A which will later become the renormalized coupling. There exist different but equivalent approaches to explain 
how the bare couplings Aq and A2 need to be chosen. Here we follow a method which, although it does not 
allow to capture the general structure behind the renormalization of ^-derivable approximations, has the advantage 
of being very similar to the approach used to renormalize the gap equation in the presence of the tadpole integral only.^^ 



The first step is to decompose the self-energy into a local and a non-local part M'^{K) = + M^^{K) with 

G(g)-Go(g)], (42) 
G{Q) G{Q - K) , (43) 



^4>' + - 
2^2 

A2,nl ,2 A^ r 



IQKA 



\Q-K\<A 



where, for convenience, we have split the bare parameter A2 as A2 = ^2,1 + A2.ni- The decomposition of the self-energy 
into local and non-local parts was used also in [l8l - l20| and it naturally arises when working in the auxiliary field 
formalism. The integral in Eq. (j43p introduces a cut-off sensitivity in the gap equation which, as we have seen, does 
not lead to a divergence of the solution but has some undesirable effects. To get rid of those, we can try to absorb 
the sensitivity by choosing 



A2,„i = AM GliQ). 
'IQKA 



(44) 



The usual consensus is that this choice is enough to absorb the sensitivity of M^^^^K) with respect to the cut-off 
because we expect that this sensitivity does not depend on the self-energy. In fact the self-energy is expected to grow 
logarithmically at large momentum. 



Let us now treat the local contribution . Using the identity: 

G{Q) = Go - {M\Q) - m^)G^{Q)G{Q) = Go - {M\Q) - m^)Gl{Q) + {M^{Q) - m' fGl{Q)G{Q) , 



we write as 



^ ^ J\Q\<K ^ 



QI<A 



{M\Q)~m'fGl{Q)G{Q). 



(45) 



(46) 



We expect the last integral not to be sensitive to large values of the cut-off. In contrast, the first integral is sensitive 
to A. Using the decomposition of M'^{Q) into a local and a non-local part in the first integral, we arrive at 



1 

A^; 



GliQ) 



<A 



A2j£ _ 1 
Aq 2 2 



Ml,{Q)Gl{Q) 



i<A 



These two different bare couplings should be seen as two different approximations to the unique bare coupling of the exact theory. 
The reason why these two bare couplings are different in a given truncation is that they will renormalize divergences originating from 
different diagrammatic topologies. 

In a forthcoming work, we shall apply the same method at finite temperature. There, we will make contact with the more general 
approaches developed in [15|| . 
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{M\Q)~m'fGl{Q)G{Q). 



(47) 



Q|<A 



We are now in a situation where we can use the same trick as in the case where the tadpole contribution is present 
only. We set 



1 



1 1 

A " 2 



Gl{Q), 



i<A 



which leads to 



Aq 



Ml,{Q)Gl{Q) 



Q|<A 



{M\Q) - m'YGi{Q)G{Q) 



(48) 



(49) 



QI<A 



However, the first integral is still potentially sensitive to the cut-ofF. We could absorb this sensitivity by adjusting 
A2,i but for this to be possible the sensitivity should be proportional to cfp' . To check this, note that M^y{K) can be 
written as 



Ml{K) ^^^^ 



A2 



Gl{Q)- \ Go{Q)Go{Q-K) 

\Q\<A "'IQI<A 
|Q--ff|<A 



Go{Q)Go{Q~K) 



IQKA 

.IQ-i^KA 



GiQ)GiQ-K) 



IQKA 

-K\<A 



(50) 



Because we expect the second term of this equation to decrease fast enough at large \K\, the sensitivity of the first 
integral in Eq. (|49p originates exclusively from the first line of Eq. (15(11) and is thus proportional to (j)^ , as needed. The 
remaining sensitivity in Eq. can thus be absorbed by adjusting A2,i such that 



GliQ) I Gl{R)- j GoiR)GoiR 

'^O ^ "'|QI<A -'|fl.|<A -'|fi|<A 



Q) 



\R-Q\<A 



= 1 



(51) 



The choice of 1 in the r.h.s. of this condition is such that A2 = Aq + O(A^), as it should be at this order of accuracy.^'^ 
Using Eq. ([5T|) and ([5(1)) in Eq. and adding the non-local contribution (H5)) with the choice of A2,ni given in 
Eq. (j44p . we finally arrive at the completely renormalized gap equation 



M^K) = + + ^ 



QI<A 



{M\Q) - m^GliQMQ) 



G{Q)G{Q-K)- Gl{Q) 

Q|<A J\Q\<A 
.\Q-K\<A 



GliQ) 



QI<A 



G{R)G{R-Q) 



fl|<A 
.l-R-QKA 



Go(i?)Go(i?-Q) 



'|-R|<A 
\R-Q\<A 



(52) 



An equation of similar form was derived using slightly different renormalization conditions in [l9j and also in |2]| . It 
was solved together with the field equation at finite temperature in Minkowski space [19i] and at zero temperature in 
Euclidean space [2l[. Here, we solve this equation for a fixed value of cf). 



We mention here that the spUtting of A2 in local and non-local parts is completely arbitrary and that the actual expression of both 
parts depend on the procedure used to determine them. Indeed, the interested reader can check that if one multiplies Eq. II47I I by Aq 
and writes its l.h.s. in a form similar to that in IIB3I I. then the inversion of the resulting operator (see | |B4|| ) results in a different form 
for <5A2.i and 5A2,nl- But this is not a problem at all, for what matters is that this different procedure gives the very same expression 
for A2 in terms of divergent integrals. 
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FIG. 11: The cut-ofi dependence of the converged value of M^(fem) obtained from the solution of the completely renormalized 
equation as compared to the solution of the toy equation with and without the inclusion of the two double integrals. 

converged with 10 ^ accuracy after 7 iterations. The inset at the left shows that M^(fcni) depends non monotonously 
on A and approaches the asymptotic value as 1/A. The inset at the right shows the approach of the converged value of M^lkm) 



in the case of the highest cutoff used. The model parameters are: 
points is [lOOA] and fe^ = 5.10"". 



1, A = 10, and 



0.2. The number of discretization 



B. Numerical solution of Eq. (|52|l and discussion 

Figure [TT] represents the numerical solution of Eq. (|52t for increasing values of A and of some approximated versions 
of this equation obtained by using a mean-field approximation and keeping or dropping the double integrals in 
Eq. (|52p. The solution of these equations are compatible with the existence of a continuum limit, which is approached 
relatively fast, as 1/A in the case of the sharp cut-off that we have chosen here. It is important to realize that the bare 
coupling redefinitions that we have obtained are precisely such that, if one would expand the solution of the original 
bare gap equation (j4T|) in powers of the renormalized coupling A, the coefficients of the expansion would converge as 
A — > c», see [2T| for a diagrammatic representation of this fact. Thus, even though the divergences of such a perturba- 
tive expansion do not always appear as divergences of the solution of the gap equation, the "correct" renormalization 
procedure needs to get rid of those as well, in order to avoid undesirable features such as those we have pointed out 
previously. This reinforces our belief that, although the two-point function is defined non-perturbatively, constructing 
a renormalization procedure which remove those divergences which appear when expanding formally the two-point 
function in powers of the coupling, is a sensible thing to do. We insist once again that this is a priori not obvious 
because the object we are finally considering is the solution of a self-consistent equation, not its perturbative expansion. 

Note that without a deeper analysis of the existence and the nature of the solutions of Eq. (|52p , we cannot really 
conclude to the existence of a continuum limit. But this is not so important, for what really matters is that there 
exists a wide range of cut-off scales, far above m and (p, where the solution exists and is almost insensitive to the 
cut-off, to within power law corrections, as it is shown in Fig. [TT] In fact from Eq. (|48|) it follows that, if one wants 
to maintain Aq > 0, one has to choose A > and A below a certain scale Ap (Landau scale) defined by 

1 1 f 1 1 1 Ap 

° " A " 2 J\Q\<A, [Q^ + rn^Y ^ A " 32^ ^ ' ^^^^ 

In the regime of interest, that is when m, (j) <^ A, Ap, we can fulfill simultaneously the requirements that the solution 
is almost insensitive to the cut-off A and that A < Ap. Taking values of A above Ap almost does not change the 
solution, although this corresponds to negative Aq. Note also that, as long as Aq > 0, we have A2 > 0. This is because 
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A2.ni > and 



A2,i — Ao 



QI<A 



GliQ) Bo{Q)-Bo{0) 



(54) 



is positive since, as we show in Appendix |d1 the perturbative bubble integral defined with Gq as 



Bo{K) 



GoiQ)Go{Q - K), 



(55) 



IQKA 

\Q-K\<A 



is a monotonously increasing function of |i^| 



Another interesting question is the following. We have considered a renormalization condition for the mass at 
(j)'^ = which somehow assumes that the system is in the symmetric phase (since (j) — is accessible). Because the 
self-energy at zero momentum should represent the curvature of the effective potential (in the exact theory at least), 
we expect it to be defined and positive for any value of 4> (at least as long as cj) does not become of the order of 
the Landau scale Ap). However this is not necessarily guaranteed because the completely renormalizcd gap equation 
involves contributions with different signs. This issue will also appear when studying the problem at finite temperature 
because we shall perform renormalization at some high temperature where we will require that the system is in its 
symmetric phase by imposing that the self-energy at zero momentum and zero field be positive. We will then have to 
check, probably numerically, that the self-energy at T = and zero momentum remains positive as we increase the 
field. Here, we propose a discussion of this issue, at T = and on a simplified version of the gap equation which allows 
for certain analytical arguments. Our main purpose is to show that M'^{0) should remain defined in some range of 
parameters, which we name the regime of interest, that is when m, </) <C A, Ap. We consider the equation 



IQKA 



QI<A 
1 
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Af 2 Q2 
1 



A'P - m2 
(g2 + m2)2 



(g2+M2)2 (g2_^^2)2 



(56) 



obtained from Eq. (j52p by performing a mean-field approximation, as in previous sections, for the zero momentum 
self-energy and by dropping the last line which contains double integrals. The solution of the completely renormalizcd 
gap equation in the mean-field approximation can be seen in Fig. Ill[ both without and with the inclusion of the 
double integrals. Interestingly, the solution obtained without the inclusion of the double integrals, that is of Eq. (|56l) . 
is closer to the solution of the original equation (1521) . For this reason we find it worthwhile to study this simplified 
version of the equation. The latter can be put in the form = /ia(M2), so that for the first and second derivatives 
one has 
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|<A {Q^ 

1 



1 



A2</)2 



+ M2)2 2 7|Q|<A (Q2+^2)2 
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A (Q2+M2)3 ' 



(g2 + IvPf 
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QI<A 



(Q2+M2)4 



(57) 
(58) 



The second derivative of hh{M^) is identical to that of gt^{M'^) which we discuss in App. [C] with the only exception 
that Ao is now replaced by A. Because A > 0, we conclude then that h'l^{M'^) vanishes only once as M"^ varies from 
to 00. It is negative for small enough M"^ , and positive for large enough M'^ . It follows that h'j^{M'^) first decreases 
from h'j^{0) = 00 and then increases towards 



h'^ioo) = -l + ^ I 

^ J\Q\<A 



(g^ 



-,2\2 



(59) 



The rest of the discussion depends on the sign of /i^(oo) and thus on the location of A with respect to the Landau 
scale Ar,: 



1. If A < Ap, /ia_(oo) < 0. Then h!f^{M'^) vanishes only once for Af2 = Af2(A) and hjy{M^) increases from 
^a(O) — to /ia(M2(A)) and then decreases towards hA{oo) = —00. This is a situation that we have already en- 
countered: the number of solutions depends on the sign of hA{M^{A)). Now, it is easily checked that h!y{m?) = A(/)2/2 
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and thus /ia(M^(A)) > 0. Then, as long as A < Ap, there are two solutions (which can be degenerate). We will be 
more particularly interested in the right most solution, which we name the "physical solution" , for it is continuously 
connected to the unique solution for A = and thus corresponds to the solution of the completely renormalized gap 
equation that we have plotted in Fig. [TTJ 

2. If A = Ap, h'jy{oo) = 0. We can draw the same conclusions as above with the only exception that, for large 
enough M^, h^iM^) does not decrease towards —oo but towards 
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2 /iQKAp 



(60) 



Then, if (j)-^ > (jy^p, with 



(61) 



the physical solution disappears. However in the regime of interest, such that in particular Ap 3> m, this situation 



corresponds to extremely large values of the field, namely > Ap/487r . In contrast, if 
solution subsists. 



< 



the physical 



3. If A > Ap, /i^(oo) > 0. There are three cases depending on the sign of the minimal value reached by h'^{M'^). 
Without trying to be exhaustive let us analyze what happens in the limit A — > oo. We have 
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(62) 



It is easily seen that h'^{AP) reaches its minimal value for — A^^, that is 



In — TT + 1 

771 



(63) 



If this quantity is positive, h'^{A4'^) > and /loo(Af^) increases from — oo to oo. There is one solution but it does 
not correspond to the physical solution. However, comparing (|63p and (I53p . we observe that this situation occurs 
only if A^^e^ > Ap, that is for extremely large values of the field. In the opposite case, h'^{AP) vanishes twice, for 
— (oo) and = M^{oo) and thus hao{Ad'^) increases from — oo to hrx>{M^ (oo)) , then decreases to ft-oo(-M^f'^(oo)) 
and finally increases again towards oo. There can be up to three solutions, the physical one corresponding to the 
"second" one. For the later to exist, we must have hoo{Mf{oo)) < 0. We have 



/Ioo(M2(oo)) 



(Aff2(oo) + A</.2)ln 



M2(oo) 



(64) 



with Mf{oo) the right-most solution of 



= -1 



327r2 



In 



M2(oo) A(/)2 



M2(oo) 



(65) 



Then, we can also write 



/ioo(M'(oo)) = m'[l + 



327r^ 
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2 327r2M2(oo) 



AAf2(oo) 
327r2 
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M2(oo) 



(66) 



Now, it is easy to check that A4f{oo) decreases as increases, from Mf{oo) ~ Ap/e when q ^ as 0^ 

approaches Ap/Ae^. Thus, in the regime of interest, M^^oo) is way larger than and (j)"^ and hoo{Mf{oo)) < 0, as 
announced. To summarize, this indicates that in the regime of interest the physical solution can be continued from 
A < Ap to arbitrary large values of A. 
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V. CONCLUDING REMARKS 



Using a particular exaniple of truncation where the two-point function has a non-trivial momentum and field 
dependence, we have revisited the renormalization of ^-derivable approximations paying particular attention to the 
question of the existence of a solution of the gap equation for arbitrarily large values of the cut-off and to the fact 
that some of the perturbative divergences which appear when expanding the solution of the gap equation in powers 
of the coupling do not appear as divergences at the level of the non-expanded solution. We have shown that it was 
nevertheless important to absorb these "perturbative" divergences in order to avoid certain inconvenient features. 

Our analysis contributes to clarify the meaning of the renormalization procedure for ^-derivable approximations 
which has been constructed in recent years. It shows in particular that the formal perturbative expansions which 
have been used sometimes to construct a renormalization scheme for ^-derivable approximations are a good guiding 
principle, although the corresponding perturbative divergences not always appear as divergences at the level of the 
solution of the gap equation. It also points to the fact that approximated renormalization schemes where (at least the 
leading) logarithmic sensitivities are not completely eliminated need to be considered with care. All these remarks go 
beyond the particular framework of ^-derivable approximations and probably apply to other resummation methods 
such as Schwinger-Dyson equations. 

Finally this work is an important step towards the inclusion of finite temperature effects. In fact, once the renormal- 
ization has been properly performed at zero temperature, the equations at finite temperature should be automatically 
renormalized. We plan to study this truncation at finite temperature in the imaginary time formalism, extending our 
approach ii i 1171 and investigate whether or not it leads to a correct order for the phase transition, as it was already 
claimed in [l9| from numerical results obtained in the real-time formalism. 
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Appendix A: General regularization for the 2PI effective action 



Consider the following generating functional 

Zii K]^ / 'D^ exp {-\ ^ ■ (Goi?)-^ ■ <f + gjntM + J-v^ + ^ifK-^} 

where J-ip = J^J{x)ip{x), ip- K ■ ip = Jyip{x)K{x,y)tp{y) and (Goi?)"^ • = Jy ip{x)G^'^ {x,y)R-\y, z)ip{z). 
For K — 0, the functional Z[J, 0] is regularized due to the presence of the regulator R and the normalization factor 
1/ J'Dip exp{— i (f ■ (GqR)^^ ■ (p}. To see this explicitly, one can consider a perturbative expansion which turns the 
evaluation of Z[J, 0] into the evaluation of Gaussian integrals. Wick theorem for Gaussian integrals generates an 
overall infinite determinant factor which is precisely canceled by the normalization factor. If follows an expansion of 
Z[J, 0] organized in terms of Feynman integrals involving the propagator GqR and thus properly regularized for an 
appropriate choice of R. In what follows, we will only consider the functional Z[J,K] in the vicinity of JsT = 0. For 
those "small" values of K, we expect the regularization at ii' = to be sufficient to regularize Z[J,K] as well. 

The 2PI effective action r[(j), G] can be seen as the double Legendre transform of Z[J, K] with respect to the sources 
J and K. Taking properly into account the presence of a regulator R and a normalization factor, one obtains 

r[cl>,G] = i0-(Goi?)-i-</)+iTr [lnG-i-ln(Goi?)-i + (Goi?)-i-G-l] +ri„t[<^,G], (A2) 

where the first two terms can be obtained by switching off all interactions in Eq. (jAip and Tint [</>, G] includes all con- 
tributions due to interactions in the form of two-particle irreducible diagrams with propagator G. The full propagator 
G[(/)] corresponding to the functional Z[J, K = 0] can be obtained from a variational principle applied to r[(/), G] that 
is 

0-^ . (A3) 
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This condition can be written equivalently as 

(A4) 



G-i[<^] = (Goi?)-i + 2^^'"* 



SG 



4>,G14>] 



Note that it is a priori not obvious that the expressions (jA2p and (jA4[) are regularized, specially because the regulator 
R does not appear explicitly in the Feynman integrals contributing to FintJ^, G]. In what follows, we address the 
question of how the regularization of Z[J, K] appears at the level of r[(/), G]. We do so in a simple situation where 
the source K{x,y) is translation invariant, allowing us to conveniently work in Fourier space. 

It is important to point out that the regularization of Z[J, K] is a priori effective when K is small enough. In terms 
of the conjugated variables this means that G is close to G. The later has a UV asymptotic behavior ^ GqR, because 
except from the momentum independent contributions originating from the tadpole integral, the quantum fluctuations 
that contribute to G are suppressed in the UV by the regulator R. It follows that in the vicinity of G, which is the 
only region where we need to consider the functional F[(/), G], the Feynman integrals contributing to rint['/>, G] are 
all regularized. It remains to discuss the trace contribution of Eq. (jA2p . Note that if one considers each of the 
contributions within brackets separately, they are all divergent. However, if we consider all terms simultaneously, we 
have 

[lnG-^-ln(Goi?)"^ + (Goi?)"^-G-l] ^ ^{UGqR)^ , (A5) 

where we introduced 11 = G^^ — {GqR)^^. Because, in the UV, 11 receives only contribution from the tadpole 
integral, it follows that the trace contribution is convergent for appropriate choices of the regulator R. 

In practice, it is convenient to consider the change of variables G — >■ GR. The gap equation then becomes 



G"M<^] =Go'+2i? 



STint 



SG 



(A6) 

4>,G[4>]R 



This is the starting point that we have considered in order to write Eq. ([T]) with the particular choice R{Q) = Q{A—\Q\). 
Remember finally that the trace contribution in Eq. (|A2p does not need regularization if all the terms are combined 
before taking the trace. One can however regularize the trace as 

Tr^ [inG^i - ln(Goi?)"^ + (Goi?)"^ • G - l] . (A7) 

This is useful for practical purposes because one can thus treat each term of the trace separately. Moreover, the gap 
equation (jA6p becomes then 



RG-^[(f>] =RGo^ +2R 



SG 



(A8) 



',G\d.]R 



In particular, the choice R = R shows that the presence of the regulator R in front of the self-energy in Eq. (jA6p is 
not crucial for the matter of regularization. 



Appendix B: Subleading behavior of the bare self-energy 



In order to determine the next term in the asymptotic expansion of M^{K) at large A, we subtract M^, in the 



form given by Eq. Q, from Eq. ([2]) with L = 0. We obtain 



G{Q)-Goo{Q) 



'^O 72 



|Q|<1 

Goo{Q)GooiQ ~ K) , 
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|Q|<1 

k\<i 



where Gao{Q) ~ ^l{Q^ + -^^i^)- Up to order contributions, one can write 
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^0 72 



Goo{Q)Goo{Q-K) + 0{ 
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(B2) 



Note that the first integral in Eq. (jB2p generates an infinite number of contributions. To resum them, we bring 
this term to the l.h.s. in the form (recall that \K\ < 1) 
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The operator that appears in the l.h.s. is invertible since 

A, 

d[L- K)- 

'\K\<1 



S(L-K)-^GUK) 



SiK^Q) + ^GUQ) 



= 5{L-Q) 



(B4) 



holds if one chooses 1/Aoo = 1/Ao + (1/2) Jj^^i^^ G'^((5). Applying the inverse operator on each side of Eq. (|B3p . 
multiplying by A^, renaming L = K/ h and using the fact that K is kept fixed, we end up with Eq. ([5]). 



Appendix C: Mean- field approximation for the mass-renormalized gap equation 

Let us here discuss the behavior as A increases of the solution of the mean-field equations (|31l) and 



1. Equation for the zero-momentum self-energy 

The mean field approximation for the zero momentum self-energy is given by Eq. Note that 50(^2) = 

-M'^+rn^ + (Ao/2)02 and thus, for A = 0, the equation = goiM'^) admits one solution only: = to^ -f (Ao/2)(/)2. 
In order to discuss the solutions of = g\{AP) as one increases A from to cx), we study the profile of the function 
gAiM"^) as A is varied. Note that its first and second derivatives with respect to are given by: 



g'AiM^) = -l- 
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(CI) 



(C2) 



It is easily checked that g'^{AP) changes sign only once. It starts infinitely negative as — > 0+, changes sign at 
some point and then remains positive, approaching as AP 00. It follows that ^^(M^) decreases from (/^(O"*") = 00 
to a certain value and then increases towards 5^(00) = —1. This means that g'j^{AI'^) vanishes and changes sign only 
once for some = M2(A). Correspondingly g^iAd"^) increases from (7a (0'*') = — 00 to g^iAl'^^K)) and then decreases 
to —00. The situation seems thus pretty similar to the one concerning the function f\{Al'^) discussed in the text: 
the existence of solutions to the equation = g\{M'^) depends on the sign of gA{Mc{-^))- However the discussion 
of the sign of gA{M^{A)) is a little bit more involved than that of /A(Af2(A)), if we want to do it for each value of 
A. Fortunately, we are only interested in what happens at higher values of A and the discussion becomes again simple. 



Consider first the equation 



A3 



^dg^A^ I 

g'idA 167r25X(A2 + M2(A))2 



A2 + A42(A) 



- 1 



(C3) 



For A large enough, the r.h.s. becomes negative. It follows that A4^{A) decreases with A and thus that AI^{A) 
converges as A — > 00. From the definition = g'j^{M^{A)) and the explicit expression (|C1I) . it follows that the only 
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possible limit for M^(A) as A — >• oo is 0. Then, multiplying the equation = g'jy{M^{K)) by Af^(A) and taking the 
limit A ^ oo, we obtain M'^{h.) In(A^) Aq^^, that is 



Af2(A) - as A-^ C50. 

° ^ ^ In A2 

We can now evaluate giy{M'^{K)) from Eq. ((33)) and study its behavior as A — > oo. We obtain 
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327r^ 
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if 7^ Aof/)^, and 



5A(M|(A))^-|^ln(lnA^), 
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if m? = Ao0^. As in the case of fA{M^{A)), the value of gA{M^{A)) at large A is driven by the presence of logarithmic 
divergences in the equation. But in contrast to /a(M^(A)), the sign of g{M^{A)) at large A now depends on the 
choice of parameters: 

1. If < Ao(/)^, gA(M^(A)) becomes negative above some Ac and the equation — g^iM"^) has no solution. This 
situation is similar to the one of equation = JaIM"^). The singular behavior — oc (Ac — A)^/^ in the 
vicinity of Ac is explained as in this later case. 

2. If > Ao0^, (/A(Afg (A)) is strictly positive for large values of A. This means that the equation = (7a(M^) has 
always two solutions for arbitrary large values of A. As we now show, the rightmost solution does not diverge as A ^ oo. 
Rather, it converges to the trivial limit — Ao(/>^. To see this let us evaluate g^iM'^) for M"^ = [m? — Ao^^) ± Ato^ 
with Am? > 0. Due to the presence of the logarithms 



gA(™' - Ao02 ± Am2) ~ T^^m^ In A^ 

327r^ 
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It follows that there exists a value of A above which gA{m? — Xo(f? — An?) > and gA{'m? — Aq^'^ + Am?) < 0. 
This means that, at least for A large enough, the rightmost solution M of the equation = gA{M'^) is such that 
— (m^ — \{)<i?)\ < Am^. Because Am^ can be taken as small as desired, this shows that the rightmost solution 
admits a continuum limit: — ^ — Xo(j?- Using this information in Eq. p3p. we obtain 
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which shows that the limit is approached very slowly. 



2. Momentum dependent equations 

We can treat Eq. (15^ along similar lines. Notice first that it can be rewritten as = gA{M'^(K); K, M^) with 
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The first and second derivatives with respect to Af2 read 
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The second derivative being negative g'jJ^M"^; K, M^) decreases strictly with from g^(0; K, M^) to g'/^{oo; K, M^) = 
— 1. If \K\ is small enough g'^{0; K, M'^) is positive and g a{M'^ ; K , M"^) increases first from gf^{0; K, M"^) up to a 
maximal value and then decreases to .gA(oo; K, A-P) = —oo. For larger values of \K\, g'{0; K, M^) could be negative 
and then g\{M'^\ K, M^) would decrease all the way from g(0; K, Af ^) to 5a(oo; K, M^) = — oo. Notice now that 

g{M'; K, M^) = ^r/.^ [b{K) ~ 5(0)] , (C12) 

where B{K) is the bubble integral, defined in the next section, where the mass in the propagators is M^. We show 
in the next section that B{K) is a monotonously increasing function of |iir|. It follows that g{M^; K, A/P) is positive. 
Then, because for > \o4>^ is defined for arbitrary large values of A, it follows that Eq. (15^ has at least one 
solution for arbitrary large values of A. Equation ([5^ is compatible with this solution having a limit as A — >■ oo. This 
limit M'^{K) obeys the equation ([55)1 . 



Appendix D: About the evaluation of integrals 
1. The bubble integral with momentum dependent self-energy 

We begin with the main integral appearing throughout this work, the bubble integral defined as 

B{K)^~ [ G{Q)G{Q-K), (Dl) 

J\Q\<A 
\Q-K\<A 

where G{Q) — 1/(Q^+M^(Q)). Working with a 4D spherical coordinate system in which the angle between the vectors 
Q and K is 6, one can do the integrals over the remaining two angles and introduce the variable P = q'^+k'^ — 2kqcos{9), 
where q = \Q\ and k = \K\ to write 



1 



A ^min((j+fc,A) 



B^^)=-^TU dqfiq) dlgil;q,k), (D2) 



with f{q) — qG{q) and g{l;q,k) — lG{l)yJ{P — (g — kY){{q + fc)^ — P). We found that an efficient and robust 
way of computing numerically this integral is achieved upon using the change of variables^** I — u{y; q) = 
^{b{q) + a{q)) + ^{b{q) — a{q)) tanh(^ sinh(y)), where b{q) = mm{q + fc. A), a{q) = \q — k\ followed hy q = v{x) = 
^ [l + tanh (-1 sinh(a;))] . After these transformations the bubble integral reads 

^W^-^JTTl/ dxJix){b{v{x))^a{v{x)))f{v{x)) dy J{y)g{u{y;v{x));v{x),k), (D3) 

where J{z) = cosij:f°l^ji^'*ij(-^)-) ■ The great advantage is that after these transformations both integrands decays with a 
double exponential rate towards both ends of the real axis, and in all practical situations it is enough to integrate 
over a relatively narrow compact interval, e.g. [—4,4], within which the integrands have a maximum. In some cases 
the adaptive numerical integration routines of the GNU Scientific Library (GSL) i24|] crashed due to the small values 
of the functions. In order to avoid this, these integrals were approximated by the integrals of the corresponding 
Chebyshev approximated functions in which the Chebyshev coefficients where calculated using a fast cosine transform 
routine of the Fastest Fourier Transform in the West (FFTW) library [2^ . We varied the number of Chebyshev poly- 
nomials between 150 and 350 and checked that the bubble integral is calculated with a relative error smaller than 10^^. 



2. The perturbative bubble integral and its monotonous behavior 

To compute the bubble integral defined with a propagator containing a momentum independent self-energy, that 
is G{Q) — 1/(Q^ + M^), we use the following geometrical picture (see also Fig. [T2|). With a 4D spherical coordinate 



First the /-integral is converted into an integral on [—1, 1] and then on the entire real axis through the Tanh-Sinh transformation, which 
is commonly used in the so-called Tanh-Sinh quadrature, see for instance [2^. |23(|. 
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FIG. 12: Region of integration for the bubble integral. The dashed circle shows a case in which 6 G [0, tt], while the dashed arc 
corresponds to a case when < < arccos a. The remaining part of the plot shows the construction used to prove that the 
bubble integral is a monotone increasing function of the momentum. 



system we have to integrate over the common region of two spheres^^ of radius A with distance k = \K\ between their 
centers. Then, the integral over two angles gives 47r, while, depending on the value of g = |Q|, the remaining angle 9 
can go from up to tt or to arccos a < tt, with a = {k"^ + — A'^)/{2kq) determined by the intersection "point" of 
the two spheres. Then, the integral of (IDip reads: 



dq 



g2 + AP 



r(g,fc,M,-l) 



dq 



A-k q^ + AP 



T{q,k,M,a) 



where 



T{q, k, M, a) 



arccos a 



d6- 



g2 + fc2 _ 2g fc cos e + Af2 ■ 
The 0- integral above can be performed analytically with the change of variable t = tan(f?/2) : 

1 



(D4) 



(D5) 



T(q,k, M,a) — -r ^ / dt- „,„ „ 

1 



where we introduced tmax — 



Aq^k"^ 
1 ~ a 



2qk^/l 



a'- 



+ (q^ + fc^ + M^) arccos a - 2e+e_ arctan | — W j ^ " 



(D6) 



1 + a 



, a 



has 



{q + kf + M2 
{q - kf + M2 



and e± — \/ {q^ kY + M"^. Note that, as special cases one 



r(g,fc,M,-l) 
T(g,fc,0,-1) 



-Max((7, k)" 



■ £+£_ 



27r 



(£+ + £-) 



2 ' 



(D7) 



Using a purely geometric argument one can prove a useful property of _Bpt(-R'), namely that it increases with \K\. 
To do so, one considers the derivative of the bubble integral with respect to K\, 



dBptjK) 
d\K\ 



2h{K) + 2h{K), 



This region is always nonvanishing if |Q|, |i<r|, |Q — i<r| < A. 
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hiK) = 



h{K) 



1 



\K\ - \Q\cos9 



'IQKA 

|Q--R-|<A 



M2 {{Q-KY + M'^Y' 



\K\ 



IQKA 



A/2 (Q - KY 



■M2 



(D8) 



and prove that it increases. Let us treat first the contribution Ii{K). The trick is to decompose the integration domain 
{\Q\ < A}r\{\Q-K\ < A} in three disjoint domains C, D and D, seeFig.HH The region D is {|Q| < K]n{\Q-K\ < 
A} n {IQjcosS > \K\}. The region D is the mirror symmetric of D with respect to the axis |(3|cos6' — \K\. The 
region C is {\Q\ < A}n{\Q- K\ < A}\{D U D). One has h = Ic + Id + Id- In region C (and also in region D), one 
has \K\ — \Q\ cos9 > 0, from which it follows that Ic > 0. In order to treat lu and Iq, for each point Q in region D, 
we introduce its mirror symmetrized Q = Q{Q). We denote by ^ = 9{Q) the corresponding angle. It is easily checked 
on Fig.IHthat for each Q e D, \K\ - |0|cos6' = -{\K\ - |Q|cos6l) > 0, \Q - K\ = \Q ~ K\ and \Q\ < \Q\. We can 
now write 



lD + In = 



\K\ ~ \Q\cose 
QeD {{Q-KY + M^Y 



1 



Q2 + M2 Q^{Q)+M^ 



> 0. 



(D9) 



Next, let us consider the contribution I2, which corresponds to the variation of Bpt{K) due to a change of the 
integration domain determined by the intersection of the two spheres (see Fig. I12p . One can see geometrically that 
when \K\ increases, the two spheres separate apart and the volume of the integration domain decreases. Therefore, 
the positive value of the integral decreases leading to the increase of Bpt{K). This argument prove the positivity of I2, 
which can be checked by a simple calculation. The support of the delta function is contained within the integration 
domain if —1 < ^ ^2qk^ ^ ^ thai is — A;| < A and q + k > A. For q,k < A the first constraint is clearly satisfied, 
so that one has 



1 



1 



327r2fc2 A2 + Af 2 



dqq- 



A2 



AP 



(DIO) 



Now, since A2 + fc2 — g2 jg positive at both ends of the g-interval and as a function of q it is strictly decreasing, it 
follows that it is positive over the whole interval. This means that I2 is positive, which completes the proof. 

We found numerically that this property of the perturbative bubble integral is inherited by the solution M (K) of 
the self-consistent gap equation containing B{K), which eventually becomes a monotone increasing function of K. 

3. The double integrals of Eq. (f52l) 



Next, we discuss the two double integrals appearing in Eq. ([52]) . After shifting the variable of integration Q — Q—R, 
one switches the order of the integrals and uses the method given in (|D 2^ . For the first double integral one obtains 



C = 



|-R|<A 



G(i?) 



'IQKA 

|Q--R|<A 



G{Q)Gl{Q-R) = -, 



1 



327r5 



drr^G{r) 



dqq G{q)—^T{q,r,m,-l) 



A— r 



dqq G{q)-—rT{q,r,m,a) 



(Dll) 



where Go{Q) = l/((9^ + m'^) and a ^ {P + q"^ - A^)/{2rq) with r = \R\ and q = \Q\. Using (|D6| and ^\ the 



derivatives can be readily done. Then, for the Q-integrals we do a Tanh-Sinh transformation as in (ID II) and evaluate 
the integrals numerically. The second double integral involves only the propagator Go and is defined as 



Go 



Go(i?) / Go(Q)G2(g-i?). 

\R\<K J\Q\<A 
\Q-R\<K 



(D12) 



It can be evaluated exactly as G, the Q-integral with a = 



— 1 can be even done analytically. 
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We shall denote by Cpt the double integral similar to C, but defined with momentum independent self-energy, that 
is with G{Q) — l/((5^ + Af^). Its divergence can be obtained by studying the divergence of its derivative with respect 
to A. Changing the order of the integration in the original form of the integral given in (|52|) one has 



da 



pt 



dK 



G{R) 



\R\<K J\Q-R\<K 

2 / 5{k^\R\)G{R) 



Gl{Q)G{Q-R)5{k-\Q\) 



R\<A 



\Q-R\<A 



Gl{Q)GiQ-R). 



(D13) 



Using in both terms the method leading to (|D 2|) one obtains 



dC, 



pt 



1 



dA 327r5 



A3 



(A2 + m2)2 Jo 



drr^G{r)T{r,A,M,a) + 



2 A3 
A2 + AP Jo 



dqq^Gi{q)T{q,K,M,P) 



(D14) 



where a = r/(2A) and (3 = g/(2A). Rescaling everything by A, one introduces 



m2/A2 and = / K"^ to 



write 



dG 



pt 



dA 



1287r5A [(1 + m2)2 



du- 



M2 



du- 



r/(";M') 



with 



f{u;M^) = uv/4~^+(u^ + M^ + l)arccos(|) 



-2 J(u2 - 1)2 + 2Af2(u2 + 1) + A/4arctan 



2-u (u + iy 



A/2 



2 + m\/ (w-l)2+M2 



(D15) 



(D16) 



We are interested in the behavior of the integrals in (jD15|) at small AP and rh'^, but because for small u one has 
f{u\lvP) ~ 7ru2/(l + AP), one can set fri^ — AP = in the second integral of (jPlsp only after subtracting from 
the integrand its leading contribution at small u, otherwise the integral would diverge in the infrared. Adding and 
subtracting this leading contribution and setting rfi^ — M"^ — whenever it is possible, one finds 



dG, 



pt 



1 



dk 



1287r5A 



du 
u 



fiu;0) 1 + 



2tt 



2ti i du- 



(^2 + to2)2 



In A 
1^ 



(D17) 



Using the symbolic manipulation program Mathematica the first integral is evaluated to tt (!), while the second integral 
can be done analytically, and for small m its leading contribution is — (l + lnm2)/2. Therefore, the divergence of 
dGpt/dA is given by (l/647r''A) ln(A/TO), which in turn means that the divergence of Cpt is (l/1287r^) ln^(A/m), 
independently on M. 



Appendix E: Numerical solution of the gap and flow equations 



The numerical algorithm used to solve the gap equation is illustrated on the completely renormalized version of 
this equation, that is Eq. (|52p. It is useful to rewrite this equation in the following form 



A2 



A 



A2 



A2 



where /c = j-ftTj and 
So 
D 



GliQ) 



1 



In 



|<A 167r^ 



A2 + 7772 



A2 



A2 



{AP{Q)-m'YGi{Q)G{Q) 



IQKA 



dq{AP{q)-m'Y<i'Gi{q)G{q) 



(El) 



(E2) 
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The integrals B{k), Co, and C are defined in (|DT]) . (|D11|) . and (ID12I) . The function M2(fc) is stored on a non- 
equidistant grid, that is at discrete values of the momentum given by k{j) = fcm + (A — k-^){^j j [N — 1)) ^ , where fc,„ is 
the modulus of the smallest momentum stored, N is the number of discretization points and j = 0, . . . , — 1. Note 
that the grid was chosen to be finer for small momenta. 

To solve the equation iteratively one starts with a constant value for M^(fc) determined by the terms on the right 
hand side of Eq. (jEl[) which contains the integrals Bq and Co (they are independent on M), then at every new 
iteration the integrals D, C, and B{K) are calculated using the stored values of M'^{K), and at the end the values of 
JVP (K) are upgraded using the right hand side of Eq. (|E1|) . One uses the one-dimensional Akima spline interpolation 
method to calculate the value of M'^{q) at arbitrary momenta, required by the integration routine used to evaluate 
the integrals. The numerical evaluation of this integrals was described in Appendix [Dl At a given value of the cut-off 
the number of discretization points used was [lOOA] or [400A/3]. 



We present below the method used to solve the flow equations pT|) and ([50)1 . Using the notations of subsection [DT] 
and working for the sake of simplicity with unsealed quantities, the flow equation (pO)) can be rewritten in the following 
form 

/•no 

dq[5{k ~ q) ^ H{k,q)]e{A- q)dAM^{q) ^ Z{k), (E3) 







where 



m = -Mt^w + ^- (e4) 



Here h{k, q) denotes the Z-integral in (|D2I) and we have also introduced 



/.A 

z{k)= I dqqG^{q)M\q)h{k,q), = dqq 



2, ^r,2, ..n ^ A ^1 ^ 3 q^ + 2AP{q) q^ + 2m^ 



{q^+NP{q)Y (q2+m2)2 



(E5) 



The equation pT|) can be written in the same form as (jE3l) . only that in both H{k,q) and Z{k) the second term is 
absent. In order to solve Eq. (|E3p numerically we store M'^{q) on a non-equidistant grid and uses the one-dimensional 
Akima spline interpolation method to calculate M'^{q) at arbitrary momentum value required by the integration 
routine used to evaluate the integrals in (|E5p . These are the only integrals remaining since upon discretization 
Eq. (jE3l) becomes 

^ [<5.,-A,i?,,]^ = Z„ (E6) 

where denotes the distance between two consecutive momenta on the non-equidistant grid. It is evident 
from (IE6p that one has to solve a system of N ordinary differential equations (ODE) of the form dMj/dA = 

Xj{MQ , . . . , M'^_i, A), where in order to know Xj the linear matrix equation X^j^o^ [^ij ~ '^iHij]xj = Zi has to 
be solved at every step of the ODE solving algorithm. To solve the ODEs we used the implicit 4th order Runge-Kutta 
method, the matrix equation was solved using LU decomposition and the number of discretization points was iV — 400 
or iV = 500 when solving Eq. jSOl) and iV = 175 when solving Eq. ([TT|) . 
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